Universal phase structure of dilute Bose gases with Rashba spin-orbit coupling 
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A Bose gas subject to a light-induced Rashba spin-orbit coupling possesses a dispersion minimum 
on a circle in momentum space; we show that kinematic constraints due to this dispersion cause 
interactions to renormahze to universal, angle-dependent values that govern the phase structure 
in the dilute-gas limit. We find that, regardless of microscopic interactions, (a) the ground state 
involves condensation at two opposite momenta, and is, in finite systems, a fragmented condensate; 
and (b) there is a nonzero-temperature instability toward the condensation of pairs of bosons. We 
discuss how our results can be reconciled with the qualitatively different mean-field phase diagram, 
which is appropriate for dense gases. 



The advent of ultracold gases has vastly increased the 
range of physically realizable many-body bosonic sys- 
tems, enabling the exploration of quantum-degenerate 
Bose gases possessing tunable interactions and band 
structure as well as internal degrees of freedom. Among 
such systems, those of particular interest involve single- 
particle Hamiltonians having degenerate ground states 
related by symmetries. Bose-Einstein condensation 
(BEG) — i.e., the macroscopic occupation of a particu- 
lar single-particle state — typically entails breaking these 
symmetries; hence the order parameter space and defects 
of such BEGs are richer than those of conventional BEGs. 
For instance, spin-1 BEGs ^ support fractionally quan- 
tized vortices, and in this sense resemble exotic fermionic 
condensates such as triplet superconductors. 

Just as these exotic defects stem from broken in- 
ternal symmetries, those of the Fulde-Ferrcll-Larkin- 
Ovchinnikov (FFLO) state, such as vortex-dislocation 
bound states [2], stem from its broken translational 
and rotational symmetries. The present work addresses 
purely bosonic analogs of the FFLO states, in which 
the degenerate single-particle ground states have distinct 
spatial wavefunctions. In particular, we consider the case 
in which the single-particle Hamiltonian possesses a dis- 
persion minimum on a circle in momentum space, so that 
BEG occurs at one or more nonzero momenta. Our work 
is motivated by a recently proposed realization of such 
a Hamiltonian, viz. a spin-i Bose gas subject to a light- 
induced Rashba spin-orbit coupling '3]. Simpler forms 
of spin-orbit coupling, having multiple discrete minima, 
have been experimentally demonstrated [4j. An alterna- 
tive approach to realizing a circular dispersion minimum 
would be to load the atoms into the excited band of an 
optical lattice; in this case, too, multiple discrete minima 
have been realized [S], and under appropriate conditions 
(e.g., "SE-even faulted" stackings of bilayer honeycomb 
lattices [6]) continuous minima are realizable. 

Spin-orbit coupled BEGs were originally addressed in 
Refs. [71 [S] as examples of unconventional condensation; 



it was argued in Ref. [5] that, for a pure Rashba cou- 
pling and isotropic interactions, a fragmented conden- 
sate should form. More recently, the case of the Rashba- 
coupled BEG was treated using mean-field theory [Hj and 
incorporating Gaussian fluctuations [T^ ; related systems 
have been studied in Refs. [TT]. In general, two phases 
have been found, depending on the spin-dependence of in- 
teractions: a time-reversal-symmetry-breaking (TRSB) 
state and a density-wave state. In the present work, 
we describe how interaction-renormalization effects qual- 
itatively change the phase diagram at low densities (see 
Fig. [T]), destabilizing the TRSB state and giving rise to a 
number-squeezed (and, in finite systems, "fragmented") 
limit of the density-wave state. These changes are due to 
the strong, emergent angle-dependence of renormalized 
interactions; such angle-dependent renormalizations are 
generic in systems whose low-energy modes occur around 
momentum-space surfaces, e.g., Fermi liquids |13| . Our 
results, while consistent with those of Ref. [F in the spe- 
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FIG. 1. Zero-temperature phase diagram as a function of the 
interaction anisotropy 7 = C2/C0 and the chemical potential 
fj,, showing the phases and transitions discussed in the main 
text. The (dashed) phase boundaries in the crossover region 
are schematic; the bold line indicates a first-order transition 
predicted by mean- field theory [9]. 
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cial case of isotropic interactions, hold more generally 
for any interactions that are repulsive in all angular- 
momentum channels. 

Our primary results are as follows. At zero temper- 
ature, we find — exploiting the properties of a quantum 
critical point introduced in Ref. [T3] — that the renor- 
malized interactions for a dilute gas universally favor a 
state in which the BEG forms at a pair of opposite mo- 
menta. For finite, weakly trapped systems, fragmented 
BEG is energetically favored over simple BEG at either 
a single momentum or a coherent momentum superpo- 
sition such as a density wave. In the thermodynamic 
limit, the fragmented BEG, though favored over a co- 
herent superposition, becomes degenerate with squeezed 
states that break translational symmetry. The resulting 
ground-state energy per particle scales unusually with 
the density n, i.e., as ri^/^; note that this scaling is the 
same as that of the "extremely anisotropic Wigner crys- 
tal" [H], which in fact approaches the fragmented state 
in the zero-density limit. At nonzero temperature, we ar- 
gue using renormalization-group (RG) methods that the 
leading instability is toward condensation of boson pairs, 
and estimate the condensation temperature. 

Model. We begin with the following effective model 
of a d-dimensional Bose gas having a circular dispersion 



minimum: 



-M+^{(|k2D|-fco)' + A:i} 



nd''fc.(^({k.})*^(ki)*t(k2)^r(k3) 



X*(ki 



*(k) 



(1) 



where ^'(ki) are Bose fields of momentum k^; k2D = 
{kx,ky); k± encodes all other momentum components; 
[/ is a possibly momentum-dependent interaction; and 
we have set h = 1. For Rashba-coupled bosons, spin- 
dependent interactions in the microscopic model imply 
momentum-dependent interactions because, for modes 
near k^, the spin is locked to the momentum. We 
shall first consider the universal properties of the gen- 
eral Hamiltonian H, and then relate these to the phases 
of the specific microscopic model considered in Ref. [3]. 
We focus primarily on the 2D case, in which fcj^ = 0, and 
then touch on the (similar) 3D case. 

We assume that energies associated with temperature 
T, chemical potential /z, system size, etc. are smaller 
than the spin-orbit coupling scale kQ/2M. Typical val- 
ues of k^^ are on the order of an optical wavelength [J], 
which is exceeded by the interparticle spacing in many 
experiments (fcQ/2M cannot be smaller than these scales 
if spin-orbit coupling is to play a significant role). 

As we are concerned with the low-energy limit, it is 
convenient to study only the degrees of freedom in a mo- 
mentum shell of thickness 2A centered on the dispersion 



minimum, giving rise to an energy scale flA = A^/2M 
intermediate between ko/2M and the low-energy scales 
fi and T. Integrating out degrees of freedom with ener- 
gies > generates effective interactions for modes with 
energies < Ha', as we discuss below, these interactions 
are further renormalized, and (for energies <^ ft a, take 
on universal values that are independent of A. A careful 
treatment of the high-energy renormalization, including 
the other Rashba bands, appears in a recent preprint |12| 
and confirms this picture. 

Quantum critical point. The model described by H has 
a quantum critical point (QGP) at /i = 0, corresponding 
to the phase transition from the empty vacuum to a BEG. 
This QGP was analyzed in Ref. |14j for fermions, but 
the analysis extends straightforwardly to bosons. Given 
that A fco, kinematics constrains the resulting form 
of the effective interaction vertices within the momen- 
tum shell (i.e., those for which all four momenta satisfy 
||ki| — fcol < A) to lie in the following channels: (i) forward 
scattering processes, which involve momentum transfer 
< A [denoted Fq^ (9) where 9 is the angle between the in- 
coming momenta]; and (ii) "Gooper-channel" processes, 
in which incoming momenta are almost equal and oppo- 
site [denoted Vqa (9) where 9 is the angle between incom- 
ing and outgoing momentum pairs (see, e.g., Ref. [13 )]. 
These channels renormalize differently: owing to the non- 
polarizability of the vacuum |16j , all renormalizations are 
due to the repeated scattering processes shown in Fig. [2^, 
which have different amplitudes in the forward-scattering 
and Gooper channels. For forward scattering, intermedi- 
ate momenta are constrained to lie in the regions shaded 
in Fig. [2j:, whereas in the Gooper channel intermediate 
momenta run over the entire circle of radius fco . 

The outcome of renormalization depends on the sign 
of the microscopic interactions. Any attractive interac- 
tions lead to an instability in the Gooper channel |14j . 
and thereby to bound states; this case is not expected 
to yield universal behavior. If, on the other hand, the 
initial interactions are all repulsive, one arrives at the 
following expressions for the renormalized interactions, 
for incoming frequencies of order ft < JIa (see Ref. (T4|): 



p (0\ _ {9) (2a) 

l + [iVm,A(0)/(27rsin0)]ln(^)' 



Fni9 = 0) 
Vn{m) 



1 + MFn^ (0) ^ko/VAin h ( 7f 

Voo{m) 

1 + MFoJm) (A:o/VaM)/; 



21^ 



, (2b) 
-,(2c) 



where fi{x) and f2{x) are scaling functions that are of 
order unity as a; — oo and approach zero as cc — )■ 0; 
and V{m) = J^"" V{9)e'"'^d9. Subscripts fl denote 
the incoming frequencies. Thus the low-energy (i.e., 
fi/r^A 0) values of all couplings are "universal," 
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FIG. 2. (a) Loop correction in the particle-particle channel, 
which governs the hierarchy of couplings at the QCP. (b) Loop 
correction in the particle-hole channel. These corrections van- 
ish at T = 0. (c) Kinematic constraints due to the dispersion 
structure: left, case of ^ 7^ 0: outgoing momenta are con- 
strained to lie in the shaded region, of area ~ A^; right, case 
of S = 0, for which the shaded region's area scales as AV^oA. 

i.e., independent of their microscopic values. Note that 
F{e = tt) = 2^ V{m). Thus, given that A/fco < 1, 
the couplings assume the following hierarchy: V^{m) ^ 
Fn{9 = tt) < Fn{9 = 0) « Fn{0 ^ 0,7r). Hence, inter- 
actions between particles at opposite momenta are neg- 
ligible compared with other interactions |17j . 

Case of T = 0. We now turn from the QCP to phases 
in its vicinity. Suppose that the system is sufficiently di- 
lute that when renormalization is cut off at a scale set 
by the chemical potential fx, the interactions are deep 
in the universal scaling regime. Then the interaction 
Hamiltonian is given by 7? ~ g, F{9 — d')ngngi, where 
ng denotes the boson density at a momentum of magni- 
tude fcp and direction 9. The hierarchy of universal cou- 
pling constants implies that H is minimized by a "frag- 
mented" state, having precisely N/2 bosons at some 6, 
and N/2 a.t 9 + tt ^8j. Fragmentation is favored owing 
to a momentum-space analog of Coulomb blockade (cf. 
Sec. 2.6 of Ref. [12]): bosons with opposite momenta do 
not interact with one another to leading order in y/ A/ko, 
whereas those at non-opposite momenta do interact. 

In more quantitative terms we can deduce the ground- 
state energy from the relation [20] /x = (n/2)Ff^{9 = 0) ~ 
(n/M)(^M/fc2)i/4^ giving 

E{N) ^ V (3) 

where A is the system area, and N± respectively denote 
the number of particles at 9 and 9 + Tr. Note that this ex- 
pression is universal, i.e., independent of the microscopic 
interaction strengths, and its unusual scaling is a conse- 
quence of the renormalization discussed above. As E{N) 
is minimized when N^ — iV_ = N/2, the ground state for 
finite TV is fragmented. Such a fragmented state can be 
understood as a density wave of wavevector fcg along the 



direction 9 with a randomly varying phase (analogous, 
e.g., to interfering independent condensates |21j). 

Fragmented states are typically unstable relative to 
simple condensates (i.e., those having a fixed phase re- 
lation) because spatial inhomogeneities tend to phase- 
lock the fragments [T^. In the present case, a phase- 
locked, coherent superposition would involve fluctuations 
of order y/N in N±, and hence cost an energy of or- 
der unity relative to the fragmented state even in the 
thermodynamic limit. Thus, a few scattering sites can- 
not overcome the tendency toward fragmentation. Sim- 
ilarly, a weak harmonic potential (i.e., of characteristic 
length much larger than the interparticle spacing) would 
not stabilize a coherent superposition relative to a frag- 
mented state, even in the thermodynamic limit, provided 
that — according to the standard prescription — the trap 
frequency a; — >■ so as to keep Nu'^ constant. This is 
because the typical matrix element between ±k due to 
the trap is of order exp(— 2fcg7V), which rapidly decreases 
as — ^ 00. 

Although a coherent superposition is disfavored in the 
thermodynamic limit, the energy cost of number fluctua- 
tions of order unity vanishes as 1/N. Thus, the ther- 
modynamic ground state (e.g., in a trap) is likely to 
be a squeezed state with small but nonvanishing phase 
variance, as opposed to the fragmented state, in which 
the phase is entirely random. This observation extends 
to translation-invariant systems, which should therefore 
exhibit spontaneously broken translational invariance in 
the thermodynamic limit. 

Implications for T — phase diagram. The dilute-limit 
phase diagram is simpler than that obtained from mean- 
field theory: it predicts that BEC occurs at two momenta 
regardless of microscopic interactions, provided these are 
repulsive. By contrast, mean-field theory predicts 
a time-reversal-symmetry breaking (TRSB) state or a 
density-wave state, depending on microscopic interac- 
tions. We now give an account of the crossover between 
universal and mean-field regimes, and estimate the min- 
imum densities required for mean-field results to apply. 

The dilute-gas results apply when, upon renormaliza- 
tion, the pertinent interactions have already achieved 
their universal forms at a length-scale shorter than the 
interparticle spacing; thus, a TRSB state is disfavored if 
-F)j(7r) < -F)t(O), regardless of whether the (larger) F^(0 ^ 
0, tt) couplings have approached their universal values. 
Note that Fa{9 = 0,7r) are related to the parameters cq 
and C2 in Ref. [3| as follows: FA(7r)/i^A(0) ~ 1 + C2/C0. 
[These relations, and similar ones for other couplings, 
can be derived as outlined following Eq. (3) in Ref. [9]. 
Provided C2 ^ Cq, all microscopic couplings are of com- 
parable magnitude.] Therefore, in terms of cq and C2, the 
TRSB state is favored only if 
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(4) 



where g is a constant of order unity. 

Note that, in addition to the TRSB phase, the Hamil- 
tonian of Ref. [S] also exhibits a regime in which a coher- 
ent superposition is lower in energy than the fragmented 
state, owing to terms of the form V'2ko''^-ko^koV'Oj which 
involve momenta of order 2fco and thus do not appear in 
H. 

These considerations lead us to the phase diagram 
shown in Fig. [T] in which there is no direct transition 
from the vacuum to the TRSB state. The transition from 
the vacuum to the density-wave state is unusual in be- 
ing a continuous transition {known to be continuous as 
the properties of the QCP are understood exactly [M]) 
at which both rotational and translational symmetry are 
broken. As a general rule (see, e.g. Refs. [5^, [21]) tran- 
sitions that break rotational and translational symme- 
try are first-order. For densities > A^, at which the 
renormalization effects discussed in the present work 
are not present, mean- field simulations show evidence of 
metastability [5]; this would suggest a first-order transi- 
tion between the density-wave and TRSB states. 

Case of T > 0. Following standard treatments of the 
dilute Bose gas [20] we assume that T ^ fi. Beyond the 
momentum scale At = 1/ V2MT, the physics is captured 
by a classical free-energy functional of the form 



5 = Jd'^k [^fi + (|k - kof] iVkP + S4 



(5) 



where S'4 denotes the set of angle- and channel-dependent 
couplings, and we have set 2M = 1. is a complex- 
field version of Brazovskii's model [22] (the relevance of 
Brazovskii's model was previously suggested in Ref. [7]). 
The initial values for the couplings in S'4 are the renor- 
malized interactions at a scale ft = T. At scales < A^, 
the vacuum is nontrivial, owing to the presence of ther- 
mal particles; hence, all couplings are renormalized by 
the particle- hole channel [Fig. |2|d]. It is convenient to 
expand F as well as V in terms of angular momenta. 
One can then implement the momentum-shell RG pro- 
cedure described in Ref. [53], by integrating out modes 
satisfying At{1 — dl) < \k — ko\ < At and rescaling 
k^ll + dl)k,TP -> [1 - {3/2)dl]'iP, and ^ ->■ n/A"^. The 
couplings transform as follows (ignoring the fiow of fi): 



dFijm) 

dl 
dVijm) 

dl 



3Fiim) 
3V, (m) 



AVi^{m) 
2(W^ 



2 (1-Mi)' 



A 



(6a) 
2 ,(6b) 



where A = 27rfco/AT. If the coupling constants at At are 
in the universal regime, one can use the fact that V <^ F 



to drop terms of order V"^. The last term in the flow 
equations drives all even V{m) (which are initially near 
zero) to negative values at some A2 — At{1 — o(At/A:o)), 
triggering a runaway growth of the even-parity V{m) cou- 
plings. Such runaway growth is associated with a pair- 
ing instability, which should in principle occur simulta- 
neously in all even-m channels. (However, as noted in 
Ref. [8j, the confining trap acts as a kinetic energy term 
of the form V^, which penalizes higher-m channels.) 

The pair-condensation temperature can be estimated 
by observing that arbitrarily weak attractive interactions 
in the Cooper channel give rise to pairs [24j whose bind- 
ing energy is A - MV'^k^. Pairing is favored for A > T. 
As T/Eo ~ (Ai/fco)' < (A/fco)2 < F ~ ln(A/Ai), one 
expects pairs to be tightly bound at length-scales compa- 
rable to I /At; at longer distances they can be treated as 
nonoverlapping. The system is thus a dilute gas of pairs, 
which condense at a temperature given implicitly [20^ by 
Tc ~ [h^n/Am) x 1/ lnln(r7,a^), where a is an effective 
dimer-dimer scattering range, which is of order A. 

Pairing would be straightforward to detect experi- 
mentally via radio-frequency spectroscopy, which should 
reveal a peak corresponding to the pair binding en- 
ergy; moreover, a pair condensate would support half- 
quantized vortices detectable via rotation. 

3D case. For this case, fcz = fc_L in Eq. (1); thus, the 
dispersion minimum is circular rather than spherical, and 
imposes the same kinematic constraints as in 2D. The 2D 
analysis thus generalizes readily; the chief difference is 
that the forward-scattering couplings in 3D renormalize 
to nonuniversal T-matrices rather than universal values, 
and the ground-state energy thus depends on microscopic 
couplings. However, Cooper-channel couplings approach 
the following universal expression as ri/f^A — ^ 0: 



Vn^l/[koM\n{n^/n)] 



(7) 



Hence, as in 2D, Y^m ^(™) ^ ^i.^ = tt) < F{e ^ tt) at 
low energies. It follows that the dilute-limit ground state 
universally preserves time-reversal symmetry. This qual- 
itative resemblance to 2D extends to the T > case, in 
which the free-energy functional — in this case, the variant 
of Brazovskii's model having two transverse dimensions 
discussed in Ref. |25) — develops a pairing instability as in 
2D. As the Cooper-channel couplings approach universal 
values more slowly, however, the conditions for the dilute 
limit to obtain are more stringent in 3D than in 2D. 
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